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ABSTRACT 

We consider the effects of inhomogeneous reionization on the distribution of galaxies at high red- 
shifts. Modulation of the formation process of the ionizing sources by large scale density modes makes 
reionization inhomogeneous and introduces a spread to the reionization times of different regions with 
the same size. After sources photo-ionize and heat these regions to a temperature > 10 4 K at different 
times, their temperatures evolve as the ionized intergalactic medium (IGM) expands. The varying 
IGM temperature makes the minimum mass of galaxies spatially non-uniform with a fluctuation am- 
plitude that increases towards small scales. These scale-dependent fluctuations modify the shape of 
the power spectrum of low-mass galaxies at high redshifts in a way that depends on the history of 
reionization. The resulting distortion of the primordial power spectrum is significantly larger than 
changes associated with uncertainties in the inflationary parameters, such as the spectral index of the 
scalar power spectrum or the running of the spectral index. Future surveys of high-redshift galaxies 
will offer a new probe of the thermal history of the IGM but might have a more limited scope in 
constraining inflation. 

Subject headings: cosmology: theory - galaxies: formation 



1. INTRODUCTION 

The reionization of the cosmic neutral hydro- 
gen, left over from the big bang, was an intrinsi- 
cally inhomogeneous p rocess (Barkan a fc Loebl 120041: 
iFurlanetto fc Loebh2005D . The basic cause of this inho- 
mogeneity is modulation by large scale structure of the 
formation process of the ionizing sources. The scatter in 
the local matter density produces corresponding scatter 
in physical observables during the epoch of reionization. 
In particular, modulation of the formation of the first 
sources by large scale density modes will introduce scat- 
ter into the redshift at which reionization is completed in 
different regions due to the scatter in local matter density 
ijWvithe fc Loebll2004|) . This implies that the redshift of 
reionization is not a single number that is universally ap- 
plicable to all regions of the universe, but rather follows 
a distribution which strongly depends on the physical 
size of the region under cons ideration l)Barkana fc Loebl 
l2004tlWvithe fc Loebll2005a|h By a fixed redshift, differ- 
ent regions of a given size expand and cool by different 
amounts since they were reheated by reionization at dif- 
ferent times. Inhomogeneous reionization leaves scars on 
the IGM in the form of a scale-dependent distribution of 
temperatures and densities. 

The fluctuations in the IGM temperature affect galaxy 
formation by increasing the minimum mass of a dark 
matter halo that would accrete gas from the photo - 
heated IGM ijHaiman et alJll996t Diikstra et all 12004 L 
Since the temperature fluctuations are scale dependent, 
the variance of the distribution of galaxies will also be 
strongly scale dependent. This will strongly influence 
the power spectrum of fluctuations in the galaxy dis- 
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tribution and alter its shape relative to the underlying 
power spectrum of the dark m atter. Previous treat- 
ments of the reh eating feedback ijBarkana fc Loeb1l200Ci 
lEfstathioulll992j) assumed that the photo-heating of the 
IGM and the corresponding suppression of galaxy forma- 
tion was uniform. In this paper we explore the implica- 
tions of the non-uniformity of this process. 

Future observations of the high redshift galaxy power 
spectrum, P(k), may be used to constrain cosmologi- 
cal parameters. Since the scale at which the dark mat- 
ter power spectrum becomes non-linear is smaller at 
high redshift, the high-redshift galaxy surveys would be 
able to probe smaller scales (i.e. modes with a higher 
wavenumber k) than their low-redshift counterparts such 
as SDSS or 2dF. The combination of the small scale data 
along with low redshift galaxy surveys and cosmic mi- 
crowave background (CMB) measurements, may signif- 
icantly improve our ability to constrain inflationary pa- 
rameters. This is particularly true for the spectral index 
n = d\ogP(k)/d\ogk, and the running of the spectral 
index a = dn/d\ogk, since the small scale observations 
increases the range of fc-modes over which these infla- 
tionary parameters can be constrained. 

The Cosmic Inflation Probe (CIP) is a proposed mis- 
sion to undertake such a h igh redshift galaxy redshift 
survey ijMelnick et al J 12005). CIP is designed to be a 
space-based large area galaxy redshift survey between 
2.5 //m and 5 /im, that would observe ~ I0 7 galaxies in 
Ha emission within the redshift interval z ~ 3 to 6.5. 
In such a survey, the galaxy power spectrum could be 
determined to an exquisite precision of < 1% over the 
wavenumber interval of k — 0.01 to 1 comoving Mpc -1 . 
Assuming linear bias these observations could potentially 
constrain the spectral index to better than An = ±0.002. 
However, we will show in this paper that the corrections 
to the galaxy power spectrum imprinted by the scale de- 
pendent fluctuations of the IGM temperature after reion- 
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ization are far greater than the existing uncertainties in 
the inflationary parameters. 

The organization of the paper is as follows. In §2 we 
discuss the manner by which large scale modes would 
introduce scatter into the redshift of reionization, and 
how this will impact galaxy formation. In §3 we present 
numerical results for the imprint of fluctuations in the 
IGM temperature on the galaxy power spectrum. Fi- 
nally, §4 summarizes our main conclusions. We adopt 
the standard AC DM cosmological mo del consistent with 
the WMAP data ijSpergel et alJl2003|) . with the parame- 
ters Q b = 0.044, n m = 0.27, n v = 0.73, n = 1, erg = 0.9 
and h — 0.72. All distances and wavenumbers through- 
out the paper are in comoving coordinates. 

2. INHOMOGENEOUS REIONIZATION 

As discussed in §1, the process of reionization was inho- 
mogeneous due to the modulation of structure formation 
by large scale density modes. The associated fluctua- 
tions in the redshift of reionization produced correspond- 
ing fluctuations in the temperature and density of the 
photo-ionized IGM. Since different regions were reion- 
ized at different times, they cool by different amounts 
and therefore obtain a distribution of temperatures at 
any redshift slice of the universe. We will assume that 
all rcionizing sources had the same spectrum and there- 
fore photo-heated the local IGM to the same tempera- 
ture. It is the distribution of reionization redshifts, not 
the scatter in the photo-heated temperatures, that is as- 
sumed to produce the distribution in IGM temperatures 
in our formulation 4 This temperature distribution pro- 
duces scatter in physical observables in the IGM, such as 
the galaxy formation rate. 

2.1. Fluctuations in the Redshift of Reionizaton 

Any region of size R in the universe has a probability 
distribution for its local matter density which will change 
the redshift at which the collapsed fraction of matter 
in the region ex ceeds some threshold and makes the 
region reionized (|Barkana fc Loefil 12004). This thresh- 
old is determined by the astrophysics of the problem, 
such as the star formation efficiency, the radiation spec- 
trum, the escape fraction of ionizing photons from galax- 
ies, a s well as the clumping f actor of the high redshift 
IGM (jFurlanetto et al.ll2004aH. However, given the red- 
shift of reionization 5 , zr, Barkana & Loeb (2004) showed 
that fluctuations around this redshift are independent of 
these astrophysical complications. The large-scale po- 
larization anisotropies of the CMB (K ogut et al.l 1200*3^ 
and th e scale of the ionized regions around quasars at 
z > 6 (IWvithe fc Loebll2004l iMesinger fe Haimanll200l 
IWvithe et al.ll2005|) indicate that the redshift of reion- 
ization is most likely between zr ~ 6 and 15. Long- 
wavelength density fluctuations contribute to the back- 
ground density and alter the rate of evolution within 

4 Other sources of temperature fluctuations, such as spatial vari- 
ations in the spectrum of the ionizing sources (resulting from vari- 
ations in the relative contributions from quasars or stars, for exam- 
ple) would only strengthen our basic conclusion that the imprint 
of inhomogeneous reionization on P(k) is far greater than the un- 
certainty associated with the inflationary initial conditions. 

5 This redshift corresponds to when an average (Sr = 0) region 
will be ionized. We assume that voids (<5r < 0) are ionized by 
nearby overdense regions at zr (so this is not the average but 
rather the completion redshift of reionization). 



a region with respect to the background. This implies 
that the time at which the collapsed fraction of the per- 
turbed region reaches a reionization threshold will be 
randomly shifted with respect to the ti me at which an 
unperturbed region reaches t he threshold l(Bardeen et alJ 
1986; Barkana fc Loebll2004t) . The collapsed fraction in 
the modulated region can be equated to the collapsed 
fraction in an unperturbed region by shifting the redshift 
of the perturbed region by 5z, 
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S c (z R + Sz) - S R 
^/2(a 2 (R min )-a 2 (R)) 



= F 



where Sr is the overdensity of the region under consid- 
eration of radius R and 5 c (z) is the collapse threshold at 
redshift z. The mass variance, namely the root-mean- 
square fluctuation in mass of regions of size R, is defined 

as 

ff2(i?) = J ^P(k)W 2 ikRl (2) 

where W{x) is the Fourier Transform of the window func- 
tion. We adopt a top-hat window function in real space 
of radius R with the Fourier Transform 



W(kR) 



kR 
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There is a minimum mass of dark matter halos into 
which the cosmic gas can colla pse, dissipate, frag- 
ment and eventually form stars ((White fc Reesl Il978t 
IRees fc Ostrikerll977j) . In fact there are two independent 
minimum masses: (i) the Jeans mass, which delineates 
the mass scale when gravitational forces exceed baryonic 
pressure support; (ii) the cooling mass, below which the 
gas cooling time is longer than the halo dynamical time 
and fragmentation into stars is not possible. Depend- 
ing on the thermal state of the IGM, either mass scale 
will determine the minimum galaxy mass. In the pre- 
reionization universe, where the IGM is extremely cold 
and the minimum mass is set by the requirement that 
the cooling time be less than the dynamic time. While 
the original reionizing sources can be formed through dis- 
sipation and fragmentation caused by H2 cooling, feed- 
back effects may have photodissociated much of the H2 
in the universe and required future dissipation and frag- 
mentation to procee d by atomic hydrogen line cooling 
l|Haiman et aLlfl997^ . The minimum virial temperature 
for atomic Hydrogen cooling is ~ 10 4 K, this consider- 
ation sets the minimum virial radius, i? m in, in equation 

Once the IGM is reionized and photo-heated to a char- 
acteristic temperature of ~ 2 x 10 4 K, the Jeans mass will 
exceed the cooling mass. In fact, we will see in the next 
section that the Jeans mass is replaced by the filtering 
mass, which incorporates the entire IGM thermal history. 
Nevertheless, we will assume that the initial minimum 
scale for ionizing sources is given by the requirement that 
the halo can through atomic Hydrogen line transitions. 

Equating the arguments in equation we can 

solve for the fluctuati on in the reionization redshift 
ijBarkana fc Loebl l2004) 
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Fig. 1. — The standard deviation of the fluctuations in the red- 
shift of reionization, <r z , as a function of region size. 



where Sq = 8 c (z)/(l + z) is approximately constant at 
high-redshift. Since 8r is a random variable drawn out 
of a Gaussian probability distribution with zero mean 
and variance a 2 {R), Sz is also a random variable. The 
relation between Sr and Sz is a linear transformation, 
and so the distribution of Sz is Gaussian, with a different 
variance. In this paper we will use Monte-Carlo sampling 
to calculate the distribution of physical quantities, such 
as reionization redshift, gas temperature and fluctuations 
in the galaxy number density. Throughout the paper 
we ensure that our mass density fluctuations are always 
physical (Sr > —I) and uncollapsed (Sr < 8 c (za)) by 
discarding regions that violate either of these conditions. 
We also assume that voids, regions with Sr < 0, will be 
photo-ionized by neighboring overdense regions; and so 
all voids will be reionized at Zr. 

Figure ^ shows the standard deviation of the fluctua- 
tions in the redshift of reionization as a function of region 
size. As the size of a region increases the scatter in reion- 
ization redshift fluctuations decreases. On scale R < 20 
Mpc the fluctuations in the redshift of reionization be- 
come substantial. 

2.2. High Redshift Galaxy Formation 

The baryons traced the temperature of the CMB with 
T cx (I + z) until they thermally de coupled from it at 
z ~ 140 (|Loeb fc Zaldarriaeal 12*004). after which they 
cooled adiabatically with T a (1 + z) 2 until the first 
galaxies formed. These galaxies produced photons above 
the Lyman limit (13.6 eV) that ionized and photo-heated 
the IGM to a temperature set by the slope of their ioniz- 
ing spectrum. If stars were the dominant source of reion- 
ization, then the IGM was heated to a temperature ~ (1- 
2) x 10 4 K (with the precise value depending on the ini- 
tial mass function of stars). After being reionized the gas 
cooled predominantly due to adiabatic cooling through 
its Hubble expansion, but at high redshift (zr > 8) the 
evolution of the gas temperature was also determined 
by Com pton cooling off the C MB and photo-ionization 
heating ijHui &: HaimanlfeOOS*) . In overdense regions, the 



gas eventually experienced adiabatic heating as it broke 
away from the Hubble flow and condensed under the lo- 
cal gravitational pull. This effect will modify our results 
only on the smallest scales. Even though filaments and 
pancakes may form at these high redshifts, the structures 
are typically not virialized. One (or two) dimensions may 
collapse early but other directions still expand and so the 
filament ary gas experiences both adiabatic heating and 
cooling ( Hui fc Gnedinl[l997j) . Numerical simulations of 
the intergalactic medium indicate that our approximate 
treatment holds for most of the gas mass (see e.g., F ig. 
1 and top left panel of Fig. 2 in lGardner et alj[2003j) . 

In our simplified treatment, we assume instantaneous 
reionization of hydrogen and ignore reionization of he- 
lium. Gradual photo-ionization of hydrogen would keep 
the IGM at a high temperature for a longer period of time 
before cooling processes begin to work. We nevertheless 
include recombinations of the ionized IGM and the heat 
input from subsequent photo-ionizations. If He I reion- 
ization occured at a similar redshift as H I, it will heat 
the IGM to a somewhat higher temperature; however, 
the temperature of the reionized IGM is uncertain by an 
even larger factor due to the unknown spectrum of the 
ionizing radiation. He II reionization is expected to heat 
the IGM at a lower redshift, z ~ 3-4 (Zheng et al. 2004; 
Theuns et al. 2002), for which the characteristic mass 
scale, M+ 6 , above which the abundance of dark matter 
halos is exponentially suppressed, is higher than the mass 
scale affected by the high IGM temperature. These cir- 
cumstances reduce the influence of He II reionization on 
the galaxy power-spectrum at low redshifts. 

Figure |3 shows the redshift evolution of the IGM tem- 
perature for several reionization redshift: zr — 6 (pur- 
ple, dot-dashed), zr = 9 (red, dotted), zr — 12 (blue, 
dashed), zr = 15 (green, long dashed) and the universal 
evolutionary track (black, solid). When reionization oc- 
curred early, the IGM temperature quickly approached a 
common evolutionary track; this effect is evident in Fig- 
ure |21 At high redshifts the timescale for photo-heating 
and Compton cooling is shorter than the age of the uni- 
verse, so the balance between these physical processes 
determines the common evolutionary track. If reioniza- 
tion occurred early enough, the z — 4 IGM temperature 
will not be sensitive to fluctuation around zr. 

Since reionization was inhomogeneous, there is a prob- 
ability distribution of possible IGM temperatures for any 
region of scale R in the universe. The actual IGM tem- 
perature of a region regulates the formation of galaxies 
in it. First, the gravitational force of the galaxy halo 
must be able to overcome the baryon pressure resist- 
ing collapse. These considerations determine the Jeans 
mass, namely the minimum halo mass into which the 
IGM gas may assemble. Second, the gas that collapses 
into a galaxy must be able to cool and dissipate in or- 
der to fragment and form stars, and so the cooling time 
of the gas must be shorter than the dynamical time of 
the halo. In the absence of molecules, this implies that 
the virial temperature of the host galaxy must be greater 
than ~ 10 4 K, so that atomic line cooling will be efficient. 

The naive Jeans mass overestimates the effects of 
baryon pressure since the matter distribution cannot in- 

6 The characteristic mass scale is implicitly defined through the 
relationship o-(M*) = 6 c (zq). 
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Fig. 2. — The redshift evolution of the IGM temperature for 
several reionization redshifts: zr = 6 (purple, dot-dashed), zr = 
9 (red, dotted), zr = 12 (blue, dashed), zr = 15 (green, long 
dashed) and the universal evolutionary track (black, solid). 



Fig. 3. — Characteristic mass scales: the Jeans mass (black, 
solid) and the filtering mass (red, dotted) at zq = 4 for several 
redshifts of reionization, zr. 



stantaneously adjust to a new thermal and hydrodynam- 
ical state of the gas ijGnedin fc HuilllQQSl) . The adjust- 
ment process occurs on the dynamical timescale of the 
halo. Incorporating the associated delay results in a 
new effective scale, the filtering scale, which replaces the 
Jeans scale and can be written in fc-space as 



1 



1 



'4(a) D(a) J a'H(a') 



da' 



da" 



(a") 3 H(a" 



(5) 

where a = 1/(1 + z) is the scale factor and D(a) is 
the growth factor of linear density perturbations. The 
baryonic power spectrum, Pb(k), which is strongly in- 
fluenced by gas dynamics on small scales is then re- 
lated to the dark matter power spectr um Pd m (k) through 
lIGnedin fc HuilToM IGnedinll2000bD 



P 6 (M) = e 4!/4W Pj m (M). 



(6) 



The thermal state of the IGM impacts structure for- 
mation through the suppression of small scale baryonic 
power and the corresponding increase in the minimum 
galaxy mass. Prior to reionization the minimum mass 
was set by the cooling time requirement. In reionized 
regions of the universe the minimum mass is now set by 
the filter mass. Figure |31 shows the Jeans mass (solid, 
black) and the filtering mass (red, dotted) at zq — 4 for 
several values of zr. 

The Jeans mass initially decreases as the redshift of 
reionization increases. The wavenumber corresponding 
to the Jeans length, 



kj(a) 



3 aH(a) 
2 c s (a) ' 



(7) 



increases for a fixed a = 1/(1 + zg) as the IGM tempera- 
ture decreases due to adiabatic cooling for small Zr. For 
large zr the IGM temperature approaches a universal 
temperature track and Jeans mass becomes a constant 
because the IGM temperature for a fixed zg is indepen- 
dent of zr. The behavior of the filtering mass is more 



complicated since the impact of the increased IGM tem- 
perature on the halo mass distribution is delayed. When 
zr is close to zg the filtering mass can be significantly 
smaller than the Jeans mass because of the afforemen- 
tioned delay. This effect can be seen in Figure|3]for small 
zr. Eventually the altered IGM temperature is capable 
of affecting the halo mass distribution and the filtering 
mass rises with increasing zr. Since the gas tempera- 
ture approaches a universal evolutionary track, the fil- 
tering mass curve becomes shallow at high zr. However, 
the filtering scale depends on the thermal history and an 
earlier reionization implies that the IGM had a high tem- 
perature (~ 10 4 K) for a longer time. This causes the 
filtering mass to monotonically increase with increasing 
zr- 

Because the temperature of the IGM is inhomoge- 
neous, its effect on galaxy formation is scale dependent 
and so it alters the shape of the power spectrum for the 
distribution of galaxies. We calculate this effect through 
Monte-Carlo simulations of the process. The fluctuations 
in the galaxy number density are defined as, 



n(R, zg) - n(z G ) 
n(z G ) 



(8) 



where n is the universal average galaxy number density, 
ignoring both fluctuations in the redshift of reionization 
and modulation of the galaxy formation process due to 
the large scale density modes. 

Using the Press-Schetcher l|Press fc Schechterl ll974T) 
model we calculate the mass function of galaxies, the 
cosmic mean of the number density of galaxies between 
mass M and M + dM formed at redshift zq , 



dn , , , . 
dM {M > ZG) 



p Sc(zg) da 2 
M V2^ct 3 dM 



exp[ 



2a 2 (M) 1 



(9) 



Then the mean number density of galaxies in the universe 
is 

n(z G )= \ dM—(M,z G ), (10) 



dM 



5 



where M m i n is the halo minimum mass corresponding to 
the filter scale, defined in equation (0). 

The ability of a halo to collapse and make a galaxy is al- 
tered by the overdensity of the large-scale region in which 
it is embbeded. The large-scale gravitational evolution 
of the region modulates the galaxy forma tion process by 
maki ng the linear barrier easier to cross l|Bardeen et alJ 
1986). Within a region of size R with mass overden- 
sity (5r, the modified mass function of halos can be cal- 
culated using the excursion set form alism, also known 
as extended Press-Sch etcher theory ijBond et alJ Il99li 
iSheth fc Leinsor] IT999) 



din 

2^(M,z G ; S R ,R) 

da 2 
X dM 6Xp[ 



p V 1 



S c (zg) - <5/ 



M V E V2^(cr 2 (M) - 
(S c (z G )-S R ) 2 
2(<7 2 (M)-a 2 (R)y 



l {R)f/ 2 
(11) 



here V L is the Lagrangian volume of the region and 
V E is its linearly evolved Eulerian volume. Non-linear 
evolution of the Eulerian volume will cause the galaxy 
number density to dramatically increase on small scales 
(R < 3 Mpc). We chose to ignore these non- linear evo- 
lutionary effects in order to clearly display the effect of 
fluctuations in the minimum halo mass, which is the fo- 
cus of this paper. 

The number density of galaxies formed in this region 
is obtained by integrating equation with the lower 
mass limit set by the region's specific mimimum mass, 



n(R, z G ) 



M„ 



dfi 

dM—(M 7 z G ;S R ,R). (12) 



The distribution function of the galaxy overdensity, de- 
fined in equation |(SJ|, can be calculated through the fol- 
lowing Monte Carlo procedure. 

1. Randomly generate the local mass overdensity (Sr) 
by sampling the probability distribution function 
as a Gaussian with zero mean and with variance 
given by <J 2 {R) as defined in equation (0. If —1 < 
Sr < S c (z G ) then proceed, otherwise resample the 
PDF. 

2. Calculate the fluctuation in the reionization red- 
shift (Sz) using equation Q. If Sz < (the re- 
gion is a void), then set Sz = 0, since neighboring 
overdense regions are expected to photo-ionize the 
voids. 

3. Assuming the region reionized a redshift zr + Sz 
and photo-heated to a temperature of 2 x 10 4 K, 
calculate the IGM temperature in the region at 
the redshift of interest zg, after taking account of 
adiabatic cooling due to Hubble expansion, Comp- 
ton cooling off the CMB and photo-heating due 
to the balance between recombinations and photo- 
ionizations of the fully-reionized hydrogen. Then 
calculate the filtering length and the correspond- 
ing minimum galaxy mass, M m i n , in the region. 

4. Use the extended Press-Schetcher formalism in 
equation ljll|l . to estimate the number density of 
galaxies with masses above M m i n in the region. 
Identify the variance of this galaxy number density 
divided by the square of the mean as 1 + u 2 G (R) . 



5. Repeat for other region sizes R. 

The with statistical moments evaluated according to 
standard definitions (see Appendix lA"|l . The variance of 
the distribution can be related to a galaxy power spec- 
trum, 



Var(S G (R)) = a 2 G (R) = 



k 2 dk 
"2^2 



■P G (k)W 2 (kR), (13) 



which can be compared to the standard dark matter 
power spectrum in order to quantify the effects of the 
IGM temperature fluctuations. 

3. RESULTS 

Based on the formalism described in the previous sec- 
tion we calculated the variance and skewness of the 
galaxy number overdensity. Figure 0] shows the mass 
variance in equation 113fl of galaxies formed at z G = 4 
for reionization redshifts of Zr — 6 (top panel) and for 
zr = 15 (bottom panel). The solid black curves are the 
Monte-Carlo simulated mass variances including the ef- 
fects of inhomogeneous reionization, while the red dashed 
curves are the linear dark matter mass variances pro- 
duced by CMBFAST 7 and adjusted for linear bias. The 
bias parameter is averaged over the mean mass function 
in equation 0, 



b = 



1 



n(z G ) 



M„ 



dfi 

dMb(M) — (M,z G ), (14) 



wher e the Eulerian linear bias is given by (IMo & White! 
1996) 

W = U»: (15) 

Oc{Z G ) 



for la rge mass halos this formula is accurate ijSheth et alJ 
|2001[) . The linear bias is calculated by expanding the 
extended Press-Schetcher mass function, equation iJHJ, 
to first order in S R . The mass averaged bias, equation 
(|14|) , can be altered to account for the fluctuations in the 
minimum halo mass due to the local mass overdensity. 
The minimum halo mass can be linearized as 



Mmin — Mr, 



M'S r 



(16) 



where M' = (dM m i n /dS R ), Mmin is the average mini- 
mum halo mass, then the mass averaged linear bias will 
have the extra term 



1 



n{z G ) 



M min 



M min + M'8 B , 



<IMb(M)^(M,z G ) 



rj -b{M mm )^—^-(M min , z G )M'S R . (17) 
n(z G ) dm 

For Zr — 6 the fluctuations in the minimum mass will 
be anti-correlated with the standard Lagrangian biasing 
effect since the minimum halo mass curve has a positive 
slope (see Fig. 0J. Physically this means that an over- 
dense region which reionized earlier than average, will 
have a larger minimum halo mass because the halo mass 
distribution has had a longer amount of time to reflect 
the photo-heated IGM gas pressure. If zr — 15 the ef- 
fects of the fluctuation in the minimum halo mass will be 

7 http://www.cmbfast.org/ 
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Fig. 4. — The mass variance in the Monte-Carlo simulation of the 
galaxy distribution (black, solid curve), Monte-Carlo simulation 
without including the fluctuations in M m ; n (blue, dotted curve) 
and the biased dark matter mass variance (red, dashed curve) at 
zq = 4. In the top panel zr = 6 and in the bottom panel zr = 15. 
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Fig. 5. — The normalized <r 2 (i?) for zr = 15 and zr = 6 (upper 
and lower solid black curves, respectively); the limiting cases of 
the best current constraints of n (red, dotted curves) and a (blue, 
dashed curves). 



weakly anti-correlated with the Eulerian linear bias. In 
the case of early reionization the IGM temperature ap- 
proaches a universal evolutionary track. This makes the 
minimum halo mass curve shallow and reduces the effect 
of inhomogeneous reionization. This feature is visible in 
r igure [H the z R = 6 mass variance is more suppressed 
than the zr — 15 mass variance on large scales where the 
linear bias is most applicable. The zr — 15 mass vari- 
ance is suppressed by approximately 10% on large scales, 
while the zr = 6 mass variance is suppressed by just less 
than a factor of 3. Additionally, the zr = 6 mass vari- 
ance drastically changes its shape near R ~ 20 Mpc. On 
this scale the fluctuations in the redshift of reionization 
begin to be substantial (see Fig. These effects would 
become more prominent, especially in the zr = 6 case, 
if we decided to consider zq = 5, since the slope of the 
M m in — zr curve would become steeper. Likewise the ef- 
fects would be more modest if we chose zq = 3. In Figure 
0|we include the output of the Monte-Carlo code ignoring 
the effect of the scale dependent fluctuations in the min- 
imum halo mass (blue, dotted curve). This comparison 
demonstrates that differences between the Monte-Carlo 
code curves and the biased dark matter mass variance 
are not dominated by numerical effects. 

The effects of radiative feedback on galaxy forma- 
tion were assumed to be restricted to the small scales 
on which individual galaxies could impact the evolu- 
tion of their neighbors. We have found that the radia- 
tive feedback produced by the process of reionization 
changes even the large scale (> 100 Mpc) amplitude 
of the galaxy mass variance (s ee Fig. St. The ioniz- 
ing sources are highly c lustered l)Furlanetto et al.i r2004b: 
IWvithe Rr. Lo7hll2005hl) . so the H II regions during reion- 
ization are much larger than the volume that an individ- 
ual galaxy could ionize. It is the large bias of the ionizing 
sources which allows the radiative feedback effects to in- 
fluence the large scale galaxy power spectrum. 

The proposed Cosmic Inflation Probe (CIP) satellite 



l)Memick et al.ll2~005D aims to measure the galaxy power 
spectrum down to k ~ 1 Mpc - in order to constrain 
inflationary parameters. Currently the best constraints 
on n and a come from a jo int analysis of WM AP data 
and SDSS Lya forest data (|Seliak et al.ll2005|) : the best 
single parameter 68.32% confidence level constraints are 
An = ±0.02 and Aa = ±0.01. Once the covariance be- 
tween n and a is inclu ded, the constrain ts considerably 
worsen (see Figure 3 in lSeliak et al.l l2005). The prospec- 
tive constraints from CIP are much finer, An = ±0.002 
l|Memick et al.ll2005D . 

In Figure we show a 2 (R), normalized to 1 at R = 
100 Mpc, for zr = 15 and zr = 6 (solid, black curves) 
and the two limiti ng cases of the best constraints from 
iSeliak et al . (2005); n = 1 ± 0.02, a = (dotted, red 
curves) and n = 1, a = ± 0.01 (dashed, blue curves). 
The effects of inhomogeneous reionization clearly dwarf 
the slight changes caused by varying the inflationary pa- 
rameters. Therefore any program which seeks to con- 
strain cosmological information from high redshift galaxy 
surveys must include the impact of inhomogeneous reion- 
ization on the shape of the galaxy power spectrum. 

The imprint of reionization on the high redshift galaxy 
power spectrum is strong because we assumed that galax- 
ies which reside in M m i n halos will be observed. If the 
survey flux limit corresponded to a much larger halo 
mass, the imprint of inhomogeneous reionization on the 
galaxy power spectrum would have been reduced 8 . How- 
ever, galaxies which are only a few times more massive 
than M m i n are still assembled through mergers of lower 
mass galaxies and their luminosity should be significantly 
affected by a depletion in the star formation rate of their 
building blocks of mass Af m j n . Additionally, supernova 
feedback may deplete the progenitors of their baryonic 
gas, and then the inhibited gas accretion may signifi- 

8 The target minimum halo mass for CIP is 3 X lO 11 A/0 
IKomatsuli2005l) . 
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cantly prevent any subsequent star formation in these 
halos (e.g. Dekel & Silk 1986). Since the effect of reion- 
ization could be larger by orders of magnitude than the 
allowed range of variations due to uncertain inflationary 
parameters, it may not be practical to confine a search 
to sufficiently massive halos where the inflationary vari- 
ations dominate. The exact effect that photo-heating 
of the IGM will have on halos with masses above M m [ n 
is difficult to estimate analytically. The analysis which 
leads to the filtering mass or similarly the Jeans mass, is 
based on linear theory. Although CIP's minimum halo 
mass is significantly above the minimum halo mass pre- 
dicted by the filtering scale, the inhomogeneous effects 
of reionization are substantially larger (see Fig. [SJ than 
the effect of the tilt and running of the primordial power 
spectrum that C1P is designed to measure. It is possible 
that inhomogeneous reionization may bias CIP's mea- 
surement of the linear galaxy power spectrum and the 
inflationary parameters. A comprehensive calculation of 
our effect requires a full hydrodynamical simulation of 
galaxy formation in a box large enough to account for 
the biasing by large scale density per turbations which 
cause inhomogeneous reionization (sec Barka na fc Loebl 
120041 for an overview of the effects of box size on simu- 
lations of reionization) . The inclusion of a homogeneous 
UV ionizing background in simulations or semi-analytic 
models of galaxy formation will not lead to the proper 
bias nor will it include the crucial effect of inhomoge- 
neous reionization which alters the shape of the galaxy 
power spectrum. 

The distortion of the galaxy power spectrum may be 
used as a new probe of inhomogeneous reionization. Even 
if a wide field survey like CIP might not have the sen- 
sitivity to observe galaxies in halos as small as M m ; n , 
the James Webb Space Telescope (JWST) 9 with ~ 1 nJ 
sensitivity most certainly will. A deep exposure (simi- 
lar in concept but deeper than the Hubble Deep Field) 
of a field with dimensions of 5 x 5 arcmin 2 will observe 
~ 10 6 galaxies between z — 3 and z = 5.5. The large 
density of galaxies (assuming source confusion is not a 
significant problem when galaxies are sliced in rcdshift) 
should eliminate the Poisson sampling errors which typ- 
ically prevent high redshift galaxy counts from placing 
strong constraints on cosmological parameters. 

4. DISCUSSION 

The process of reionization is intrinsically inhomoge- 
neous due to the modulation introduced by large scale 
structure on the formation of ionizing sources. Regions 
which are overdense form structure earlier than average 
and correspondingly reionize earlier than average. As 
the regions are photo-ionized they get photo-heated to 
a characteristic temperature of ~ 2 x 10 4 K, and sub- 
sequently cool. Regions that reionized earlier than aver- 

9 http://www.jwst.nasa.gov/ 



age experienced suppressed galaxy formation for a longer 
amount of time. Consequently, the minimum mass of 
dark matter halos that can accomdate star formation, 
M m i ni is higher in such regions. This feedback of reion- 
ization on the minimum halo mass alters the galaxy 
power spectrum. Figure summarizes the effects of in- 
homogeneous reionization on the distribution of galaxies 
at high redshifts. 

The filtering scale differs from the Jeans scale since 
it involves an average over the thermal history. Imme- 
diately following the photo-heating of the IGM, the fil- 
tering scale lags behind the Jeans scale. The delay is 
displayed in Figure where the filtering mass rises ini- 
tially rapidly and later more gradually as the redshift of 
reionization increases, while the Jeans mass declines at 
first and then approaches a constant. If reionization oc- 
curred at zr = 6, the overdense regions which reionize 
earlier would have a larger M m i n than average, reducing 
the number density of galaxies. This is anti-correlated 
with the well-known biasing effect of large scale den- 
sity modes which produces more galaxies in overdense 
regions. If reionization occured at Zr — 15, the marginal 
increase in the minimum mass for overdense regions is 
more modest since the photo-heated IGM quickly ap- 
proached a universal evolutionary track 

The new effect we calculated is scale dependent, since 
the fluctuations in the redshift of reionization are larger 
on small scales (see Fig. QJ. This consequently implies 
that the proper interpretation of the high redshift galaxy 
power spectrum will require a detailed understanding of 
reionization. Without adequate knowledge of the IGM 
temperature history, it may not be possible to constrain 
the inflationary parameters (through the spectral index 
and the running of the spectral index of scalar perturba- 
tions) to a high precision. Conversely, the shape of the 
galaxy power spectrum at high redshifts contains crucial 
information about the scale dependence of the thermal 
history of the IGM during and after the epoch of reion- 
ization. Previous attempts to measure the IGM tem- 
perature through its imprint on the small scale p ower 
spectrum of the Lya forest (see iZaldarriagal 2002 , and 
refere nces within) are not easy to implement l)Lai et all 
2005). Since the process of galaxy formation is inherently 
non-linear, the IGM temperature flucuations are ampli- 
fied and become easier to detect. A deep exposure by 
JWST or a wide field galaxy survey like CIP may prove 
to be effective new probes of the epoch of reionization. 
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APPENDIX 

MOMENTS ESTIMATORS 

In this work we Monte-Carlo several different physical quantities (mass overdensity, reionization redshift fluctuation, 
filtering length and galaxy overdensity) in order to determine the quantities' probability distribution function (PDF). 
The PDFs can be completed described by the infinite set of their statistical moments. In this work we will only 



measure the lowest three moments. The first moment, the mean, is defined as 

Mean(xi) = x = — ^2 Xi > C^- 1 ) 

i 

the variance is defined as 

Var{ Xi ) = a 2 = j±- - xf. (A2) 

i 
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